Defect modes in one-dimensional photonic lattices 
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Linear defect modes in one-dimensional photonic lattices are studied theoretically. For negative 
(repulsive) defects, various localized defect modes are found. The strongest confinement of the 
defect modes appear when the lattice intensity at the defect site is non-zero rather than zero. When 
launched at small angles into such a defect site of the lattice, a Gaussian beam can be trapped and 
undergo snake oscillations under appropriate conditions. OCIS numbers: 190.5530, 230.3990. 



Light propagation in periodic photonic lattices is under 
intensive study these days due to their novel physics and 
light-routing applications [1,2]. Most of these studies fo- 
cused on nonlinear light behaviors in uniformly periodic 
lattices [3-10]. A natural question arises: how does light 
propagate if the photonic lattice has a local defect? In 
photonic crystals, this question has been much analyzed 
[11]. While nonuniform arrays of "fabricated" waveg- 
uides with structured defects were used in previous stud- 
ies [12-14], the issue of defect modes in optically-induced 
photonic lattices has not yet received much attention. 
Since photonic lattices differ from photonic crystals on 
many aspects (for instance, the refractive-index variation 
in an induced photonic lattice is typically several orders 
of magnitude smaller than that in a photonic crystal), 
one wonders if photonic lattices with a local defect can 
also support defect modes. 

In this Letter, we theoretically analyze linear defect 
modes in one-dimensional photonic lattices with a local 
negative defect as induced in a biased photorefractive 
crystal. In such a defect, the lattice intensity is lower 
than that at nearby sites (akin to an "air defect" in pho- 
tonic crystals [11]), thus light has a tendency to escape 
from the defect to nearby sites. However, we found that 
localized defect modes do exist due to repeated Bragg 
reflections. More interestingly, strongly confined defect 
modes appear when the lattice intensity at the defect 
site is non-zero rather than zero. As the lattice poten- 
tial increases (by raising the bias field), defect modes 
move from lower bandgaps to higher ones. If a Gaussian 
beam is launched at small angles into the defect, it can 
be trapped and undergo robust snake oscillations inside 
the defect site without much radiation. 

The physical situation we consider here is that an or- 
dinarily polarized lattice beam with a single-site neg- 
ative defect is launched into a photorefractive crystal. 
This defected lattice beam is assumed to be uniform 
along the direction of propagation. Meanwhile, an extra- 
ordinarily polarized probe beam with a very low intensity 
is launched into the defect site, propagating collinearly 
with the lattice beam. The non-dimensionalized model 
equation for the probe beam is [5] 
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Here U is the slowly-varying amplitude of the probe 
beam, z is the propagation distance (in units of 
2k\D 2 1 / 7r 2 ), x is the transverse distance (in units oiD/ir), 
Eo is the applied DC field [in units of tt 2 l(k^n\D 2 r '33)], 
II = Iq cos 2 x {1 + efo(x)} is the intensity function of 
the photorefractive lattice (normalized by the dark irra- 
diance of the crystal Id), Io is the peak intensity of the 
otherwise uniform photonic lattice (i.e., far away from the 
defect site), fo{x) is a localized function describing the 
shape of the defect, e controls the strength of the defect, 
D is the lattice spacing, ho = 2tt/\o is the wavenumber 
(Ao is the wavelength), k\ = kon e , n e is the unperturbed 
refractive index, and r33 is the electro-optic coefficient 
of the crystal. In this letter, we assume that the de- 
fect is restricted to a single lattice site at x = 0. Thus, 
we take Jd{x) = exp(— x 8 /12S). Other choices of defect 
functions fp give similar results. When e < 0, the light 
intensity II at the defect site is lower than that at the 
surrounding sites. This is called a negative (repulsive) 
defect where light tends to escape to nearby lattice sites. 
For e = — 1 and —0.5, the corresponding lattice intensity 
profiles are displayed in Figs, la and 3b respectively. 
In the former case, there is no light at the defect site, 
while in the latter case, there is still light at the defect 
site but with a halfway reduced intensity. These lattices 
with structured defects might be generated experimen- 
tally by optical induction. Consistent with our previous 
experiments [7], we choose parameters as follows: the 
lattice intensity Io = 3Id, lattice spacing D = 20/im, 
Ao = 0.5/im, n e = 2.3, and r33 = 280pm/V. Then one 
x unit corresponds to 6.4/im, one z unit corresponds to 
2.3 mm, and one Eo unit corresponds to 20 V/mm in 
physical units. 

For a negative defect, a surprising feature is the pos- 
sible existence of localized defect modes due to repeated 
Bragg reflections. The existence of such modes will af- 
fect light propagation in a profound way. We seek such 
modes in the form U(x,z) = e~ lfiz u(x), where function 
u(x) is localized in x, and \i is a propagation constant. 
Our numerical method is to expand the solution u(x) 
into discrete Fourier series, then converting the linear 
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u(x) equation into an eigenvalue problem with \i as the 
eigenvalue. First, we consider the defect with e = — 1, 
where the lattice intensity at the defect is zero (see Fig. 
la). For this defect, we have found defect modes at vari- 
ous values of E$. The results are shown in Fig. lb. It is 
seen that at low values of Eo (low potential), two defect 
modes appear in the first and second bandgaps. The one 
in the first bandgap is symmetric in x, while the one in 
the second bandgap is anti-symmetric in x. Both types of 
defect modes are moderately confined. Examples of such 
modes at Eo = 1.2 and 3.5 are displayed in Fig. l(c, 
d) respectively. However, these defect modes disappear 
when Eo increases above certain threshold values. In par- 
ticular, the symmetric branch in the first bandgap disap- 
pears when Eo > 2.8, while the anti-symmetric branch in 
the second bandgap disappears when Eo > 7.5. On the 
other hand, before the antisymmetric branch disappears, 
another symmetric branch of defect modes appears in- 
side the same (second) bandgap. This new branch exists 
when 5.3 < Eo < 10.3, and it is generally more localized 
than the previous two branches. This can be seen in Fig. 
le where this symmetric defect mode at Eo = 7.5 is il- 
lustrated. Compared to Fig. l(c, d), this new mode is 
much more confined. 

Three general features in Fig. 1(b) should be noted. 
First, for any positive Eo value, at least one defect mode 
can be found. Second, each branch of defect modes dis- 
appears as Eo increases to above a certain threshold. 
Thirdly, as Eo increases, defect modes disappear from 
lower bandgaps and appear in higher bandgaps. In other 
words, defect modes move from lower bandgaps to higher 
ones as Eo increases. 

The existence of these defect modes as well as their 
profile and symmetry properties have a profound effect 
on linear light propagation in the underlying defected 
photonic lattices. If the input probe beam takes the pro- 
file of a defect mode, then it will propagate stationarily 
and not diffract at all. This is seen in Fig. 2(b), where 
the numerical evolution of an initial defect mode (with 
e = — 1 and Eo = 7.5) is displayed (the corresponding 
lattice field is shown in Fig. 2a). For a Gaussian in- 
put beam (as is customary in experimental conditions), 
the evolution will critically depend on whether a defect 
mode resembling the input Gaussian beam exists under 
the same physical conditions. To demonstrate, we take 
an initial Gaussian beam as U(x,0) = e~^ x which re- 
sembles the central hump of the defect mode in Fig. le, 
and simulate its evolution under various Eo values by 
pseudo-spectral methods. The lattice intensity field is 
the same as that in Fig. 2a (where e = —1). We found 
that at small values of Eo, the Gaussian beam strongly 
diffracts and quickly becomes invisible. Similar behav- 
ior persists as Eq increases (see Fig. 2c) until it reaches 
a value about 7.5, when a large portion of the initial 
beam's energy is trapped inside the defect site and prop- 
agates stationarily (see Fig. 2d). As Eo increases beyond 



7.5, however, strong diffraction of the probe is seen again 
(see Fig. 2e). These results indicate that the light trap- 
ping in Fig. 2d could not be attributed to either the 
simple guidance due to increased lattice potential or the 
nonlinear self-action of the probe beam itself. Rather it 
must be attributed to the repeated Bragg reflections in- 
side the photonic lattice under certain phase-matching 
conditions, as the Gaussion beam matches the localized 
mode of the defect. This bears strong resemblance to 
localized modes in photonic crystal fibers. 

For various applications, it is often desirable to keep 
the defect modes as locally confined as possible. The de- 
fect considered above with e = — 1 (see Figs, la and 
2a) is certainly simple and intuitive, but does it give 
the most strongly confined defect modes? To answer 
this question, we fix the value of Eo and allow the de- 
fect parameter e to vary from —1 to 0, then determine 
at what e values the most localized defect modes arise. 
With fixed Eo = 6, we have obtained the defect modes 
versus e and plotted the results in Fig. 3. Fig. 3(a) 
reveals that at small negative values of e, a single de- 
fect mode bifurcates from an edge of a Bloch band inside 
each bandgap. As e decreases, the defect mode in the 
first bandgap disappears (at e = —0.81), while the one 
in the second bandgap persists. The defect-mode branch 
in the first bandgap is more localized than the one in 
the second bandgap in general. Thus we focus on this 
branch in the first bandgap below. When |e| is small, the 
defect eigenvalue is rather close to the left Bloch band, 
thus the defect mode is rather weakly confined (see Fig. 
3c). As |e| increases, the mode becomes more confined. 
As e approaches —0.81, the defect eigenvalue approaches 
the right Bloch band, and the defect mode becomes less 
confined again (see Fig. 3e). Surprisingly, we found that 
the strongly confined defect mode occurs when e w —0.5. 
This defect mode and the corresponding lattice intensity 
field are shown in Fig. 3(d, b) respectively. These find- 
ings are rather interesting, as they show that the most 
localized defect mode arises when the lattice intensity at 
the defect site is non-zero rather than zero. Such results 
may have important implications for applications of de- 
fect modes in photonic lattices. 

We have further studied the evolution of a Gaussian 
input beam launched at small angles into a photonic lat- 
tice with Eo = 6 and e = —0.5. For this purpose, we take 
the initial condition as U(x,0) = e~^ x + lkx ? where this 
Gaussian intensity profile resembles the central hump in 
the defect mode of Fig. 3d, and the phase gradient k is 
proportional to the launch angle of the Gaussian beam. 
At zero launch angle (k = 0), a vast majority of the input- 
beam's energy is trapped inside the defect and propa- 
gates stationarily (see Fig. 4b). When compared to Fig. 
2, we see that the confinement of the probe beam by 
the present defect (shown in Fig. 4a) is more efficient, 
mainly because the defect mode admitted under these 
conditions is more localized (see Fig. 3d). Next we take 
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k = 1, which corresponds to a launch angle of 0.58° with 
physical parameters listed earlier. In this case, most of 
the light is still trapped inside the defect site. However, 
the trapped light undergoes robust snake-like oscillations 
as it propagates through the defect (see Fig. 4c). The 
ability of a negative defect to trap oscillating light beams 
is a remarkable feature that merits further investigation. 



In summary, we have analyzed linear defect modes in 
one-dimensional photonic lattices with negative local de- 
fects. These results are expected to pave the way for ex- 
perimental observations of such localized modes as well 
as for the study of nonlinear defect modes. 
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FIG. 1. (a) The lattice intensity profile with Iq = 3 
and e = —1; (b) the applied dc field parameter 
Eo versus the defect eigenvalues /x; the shaded re- 
gions are Bloch bands; (c, d, e) three defect modes 
at (E ,n) = (1.2, 1.604), (3.5, 5.812), (7.5, 7.997) which are 
marked by circles in (b) respectively. The shaded stripes in- 
dicate the locations of the lattice's peak intensities. 
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FIG. 2. (a) The lattice intensity field with J = 3 and 
e = —1; (b) evolution of an exact defect mode (shown in Fig. 
le) at Eo = 7.5; (c, d, e) evolutions of a Gaussian beam at 
three Eq values 5, 7.5 and 10 respectively. 
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FIG. 3. (a) The defect strength e versus the defect eigenval- 
ues /x; (b) intensity profile Il(x) of the photonic lattice with 
e = —0.5; (c, d, e) three defect modes of the first bandgap 
with (e, n) as marked by circles in (a) respectively. 
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FIG. 4. Evolution of a Gaussian beam launched at zero (b) 
and non-zero (c) angles into the defect site of a photonic lat- 
tice (a). Intensity fields are shown. Here io = 3,i?o = 6 and 
e = —0.5 in Eq. (1). The initial phase gradient in (c) is k = 1. 
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